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O ■ 1 Introduction 

As hadrons can be viewed as excitations of tlie QCD vacuum, tlie vacuum structure that features 
I varieties of quark and gluon condensates plays a crucial role in the hadron dynamics [1, 2, 3, 4, 5]. 

, The chiral or quark condensate < qq > serves as the order parameter for the spontaneous chiral 

CN ■ symmetry breaking and reflects the nature of the QCD vacuum. For a chirally invariant vac- 

1^ _ uum, the chiral condensate vanishes. The non-vanishing chiral condensate plays a crucial role 

, in generating the dynamical chiral symmetry breaking with the mass acquisition of constituent 

■ quarks. In free space, < qq >o= —(225 it 25 MeV)^ according to the non-perturbative QCD ap- 

proaches, for instance, see Refs. [4, 6, 7], whereas the perturbative vacuum is chirally symmetric 
for massless quarks. In the hot and/or dense medium, the partial restoration of chiral symmetry 
with reduced in-medium < qq > is expected. The hadronic and electromagnetic signals for 
the partial restoration of chiral symmetry have been suggested in nucleus- nucleus reactions, for 
instance, see [3, 8, 9]. 

Due to the difficulties of the non-perturbative QCD, the in-medium chiral condensate may 
be derived alternatively using the QCD sum rules [10], and/or effective QCD models [1, 3, 11, 
12, 13, 14]. The Nambu-Jona-Lasinio models and the linear a models are among the well- 
known examples of the effective QCD models. Recently, the in-medium chiral condensate was 
also exploited in some new attempts [15, 16, 17]. A reduced in-medium chiral condensate 
was obtained in these models and approaches. However, an appropriate description of the 
nucleonic medium in terms of the quark degree of freedom is still very difficult. Interestingly, 
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a model-independent approach to evaluate the in-medium chiral condensate was developed by 
applying the Hellmann-Feynman (HF) theorem [11]. On the hadron level, the in- medium chiral 
condensate can be obtained from the derivative of the nuclear matter energy density with respect 
to the current quark mass. In fact, the in-medium chiral condensate has been evaluated using 
the nuclear matter energy density obtained from various many-body theories. For instance, 
the Dirac-Brueckner (DB) approach [18, 19], the relativistic mean field (RMF) models [20, 21, 
22, 23], the quark meson coupling models [24, 25], and the effective field theory [26, 27], have 
all been used in studying the in-medium chiral condensate. Surprisingly, most of these studies 
have shown various hindrances towards the chiral symmetry restoration at high densities when 
nucleon-nucleon interactions are taken into account even in the modified quark-meson coupling 
model [28]. 

In principle, the non-vanishing chiral condensate arises from the vacuum of the non-perturbative 
QCD. In particular, the nucleonic vacuum can be considered as one of hierarchic constructions 
of the non-perturbative QCD vacuum. At finite density, the vacuum is modified by the medium. 
In dealing with the chiral condensate in nuclear medium, a proper treatment of the nucleonic 
vacuum is thus very important. However, so far very little attention has been paid to the vacuum 
corrections in evaluating the in-medium chiral condensate because usually only effective models 
are used. In this work, we explore the in-medium chiral condensate with the vacuum corrections 
using the Lagrangian constructed by us recently [29, 30] adopting the Brown-Rho (BR) scaling 
law for the partial restoration of chiral symmetry in nuclear matter [31, 32]. 

The paper is organized as follows. In Section 2, we rcnormalize our model with the Brown- 
Rho scaling in the Hartree approximation to attain the finite contribution from the Dirac sea. In 
section 3, we present the formalism for the in-medium chiral condensate using the HF theorem. 
Results on the Equation of State (EOS) and in-medium chiral condensate are discussed in Section 
4. The summary is given in Section 5. 

2 Relativistic Hartree model with BR scaling 

In the present work, the model Lagrangian with density-dependent couplings and meson masses 
is written as [30] 

L = i^iij^d^^ -M + g*a- g^^^u^ - 5;7mT36^] V + \ [d^od'^a - mfa^) 

-\f^,F^^ + \m*^u^u^^ - \b,.B>^^ + \mfh^X - \a,,A^^ + Cct, (1) 

where Cct is the counter term used for renormalization, ip,a,LO, and bo are the fields of the 
nucleon, scalar, vector, and isovector-vector mesons, with their masses M,m'^,m1^, and m*, 
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respectively. The meson coupling constants and masses with asterisks denote their density 
dependence, given by the BR scaling [29, 30, 33]. The density dependence of parameters is 
described by the scaling functions that are the ratios of the in-medium parameters to those in 
the free space. We take the scaling function for meson masses as given in Ref. [29] 

^{p) = l-yp/po, (2) 

where the linear dependence on density is similar to the Nambu scaling [32] but with a much 
smaller coefficient y 0.1 suggested by recent experiments [34, 35]. Here, the scaling for the 
nucleon is not considered since the nucleon mass already drops in the medium due to the coupling 
to the a meson. The scaling functions for the coupling constants of scalar and vector mesons 
read 

^aKP) = 7- ^> ^pKP) = r— 1—, ^oj{p) = — 1—, (3) 

^ + xp/po 1 + ypp/po i + yojP/po 

which are necessary for an appropriate description of nuclear matter properties [29] . 

The model renormalization is treated in the Hartree approximation. The Lagrangian for the 
counter terms reads 

£cT = aa + ^f5a^ + ^ja^ + ^Aa^ (4) 

where the definition and meaning of the coefficients in Cqt follow those given in Refs. [36, 37] 
but with the replaced by g*. These counter terms are used to renormalize the nucleon scalar 
self-energy, while no renormalization procedure is needed for the nucleon vector self-energy due 
to the nucleon current conservation. We use the procedure developed by Chin in Ref. [37] 
to renormalize the scalar self-energy and the energy density. In the following, we just stress 
the renormalization for the rearrangement term. The latter is essential for the thermodynamic 
consistency in deriving the pressure. In the Hartree approximation, the expectation value of the 
rearrangement term T,q is given by 

^° dp ^~ ^^"^ dp dp "^'^'^ dp^^ ' 

T , dgl da I dB ^ 1 d^ n 1 d\ a 

i=p,n ^ ^ i=p,n ^ ' 

Here = g^/mlj, Cp = g*/m*, 5 = {pp—pn)/p is the isospin asymmetry, ps is the scalar density, 
and G{k) is the nucleon propagator. After the renormalization, the T,'-'' can be expressed as 



and 



with 



i=p,n 



(8) 
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where the scalar nucleon self-energy S5 = M - M* with M* = M - g*a = M - g*^ps/m*'^. 
The energy density and pressure read, respectively, 

£ = Icy + Icy 5^ + i(7^(M* - M)^ + ^ E* + AS (9) 

i=p,n V / 

P = \cy ^\cyi' - \cl{M^ -Mf -Y.^ p 

+- y / A — - AS, (10) 

where C^- = m^/g*, and is the finite vacuum correction [36, 37] 

= (M*^ln^ + m3S5 - ^M^Sl + f MS| - . (11) 

3 In-medium chiral condensate 

The HF theorem transmits a parameter dependence of the Hamiltonian to that of the system 
energy. For the QCD Hamiltonian, one may divide it into a chirally invariant part and a part 
depending on the current quark mass, namely Hqcd = Wo + fnqqq. By taking mq as the very 
parameter, the HF theorem reads 

Since the vacuum expectation value (VEV) < Q\Hqcd\^ > is subtracted to obtain the total 
energy of nuclear matter, the in-medium chiral condensate < qq >p is given as [11] 

dS 

< qq >p=< qq >o (13) 

Here, S is given in Eq.(9), and can be written more generally as 

S = {M + E/A)p, (14) 

where E/A is the binding energy per nucleon. Both hadron masses and meson-nucleon coupling 
constants depend on the current quark mass. However, based on an analysis in the linear a 
model, the derivative of the coupling constant with respect to the current quark mass is smaller 
than that of the hadron masses by a factor ctn/M [19], with the pion-nucleon sigma term 
gm ~ 45 MeV. Thus, we neglect the dependence of the coupling constants on the current quark 
mass. The in-medium chiral condensate is then given by 

/ dM d(E/A) dM ^ d(E/A)dm*\ 
<qq>p = <^-^>0+( ^ + ^^E^^)P- (15) 
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The current quaxk mass derivative of the nucleon mass is expressed in terms of the as [11, 18] 

dM _ aN_ 
druq ruq 

For the vector (a; and p) mesons, one may resort to the constituent quark model and give the 

vector meson sigma term using the following relation 

rriq duly _ rriq dM _ on_ 
my drUq M drUq M' 

This approximation has been widely used in the literatures [11, 18, 19, 28]. The determination 
from the ttN scattering gives (jjv = 45 ± 7 MeV [38]. This central value can be reached in the 
constituent quark model when the quark confinement condition is imposed [13]. Recently, in 
Ref. [27] the two-pion exchange processes including the intermediate A isobar excitations instead 
of the surrogated cr-meson exchange were considered, and a hindered tendency toward the chiral 
symmetry restoration at high densities was still shown. In these chiral effective field theories [26, 
27] , the short-range Nucleon-Nucleon (NN) interactions only play a minor role on the change of 
the in-medium quark condensate, while the long- and intermediate-range NN interactions from 
one- and two-pion exchanges dominate the hindrance. This observation is instructive, though 
different from that in the one-boson-exchange approaches [18, 19, 20, 21, 22, 23, 24, 28]. However, 
the strict renormalization for the low-energy physics and contact interactions in these chiral 
effective field theories is nevertheless imperative to further judge whether the additional degrees 
of freedom are necessary for describing the short-range interactions. Though the progress of 
renormalization has been made recently [39, 40, 41], this is still quite challenging. In these chiral 
effective theories that were built upon the chiral pertubative theory, the vacuum is alienated 
in two folds: First, the pion as the Goldstonc boson that resides in the coset space is not 
accompanied by other boson degrees of freedom that carry the chirally invariant vacuum in 
the conserved subgroup space. Second, in a non-relativistic treatment the nucleonic vacuum is 
out of touch. In fact, suitable renormalization scheme of the baryonic vacuum allows one to 
include consistently the vector mesons as explicit degrees of freedom in the Lorcntz-invariant 
baryon chiral perturbative theory [42]. Since the chiral symmetry restoration is expected to 
occur at the short-range region with a chirally invariant vacuum, the present approach is of 
close relevance to dealing with the in-medium chiral condensate regarding the importance of 
the vacuum correction and the important role of short-range w-meson exchange. As we focus 
on the vacuum correction, the relation (17), though phenomenological, can provide appropriate 
and referential information. 

The situation with the scalar a meson is more involved since it has a broad mass spectrum 
starting roughly above 2m7r, as a result of the equivalence of a variety of two-pion exchanges. 
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As known in the linear a model of Gell-Mann and Levy [43], the a meson is regarded as the 
chiral partner of the Goldstone boson (tt meson). In this model the chiral symmetry and its 
spontaneous breaking are well manifested. Although one can have some other choices to scale 
the sigma term for the a meson [11, 18, 19, 28], we stress here that the symmetry of chirality 
is important and thus in our model that is designed to respect the chiral symmetry restoration 
at high densities we persist in using the sigma term for the a meson obtained from the linear a 
model, and at tree level it is [12, 19, 44] 

ml druq 2 M ■ ^ ' 

The use of this relation is consistent with present relativistic Hartree approach which is also 
at tree level, though there can be high-order corrections to this relation. Here, the cr-meson 
mass is taken to be 590 MeV. It was obtained by fitting ground-state properties of many nuclei 
with the same scaling functions used in our previous work [30]. Using relations (2), (16), (17), 
(18) and the Gell-Mann-Oakes-Renner relation m^f^ = —niq < qq >o [6], the in- medium chiral 
condensate is eventually given by 

<gg>P _ i_ / d{E/A) d{E/A) 3m , p ^p5^\ 

with the TT-meson decay parameter = 93 MeV and the 7r-meson mass m.,^ = 138 MeV. The first 
term in the bracket is the leading-order model-independent contribution. In this work, we focus 
on the chiral condensate in symmetric nuclear matter where the p meson does not contribute. We 
also do not consider the (5, tt and ry mesons, as they do not contribute to the symmetric nuclear 
matter either in the Hartree approximation. In fact, within the Dirac-Brueckner approach, their 
contributions to the in-medium condensate are about two orders of magnitude smaller than 
those of (7 and oj mesons [18, 19]. Note that the role of the tt meson in these relativistic one- 
boson-exchange theories is different from that in the effective field theories [26, 27], where the 
perturbative expansion with the unique tt meson degrees of freedom results in the multi-pion 
exchanges [39]. 

4 Results and discussions 

Before presenting our results on the in-medium chiral condensate, we discuss more about the 
determinations of model parameters. In our previous work [29, 30], besides the saturation 
properties of nuclear matter, namely, the incompressibility k = 230 MeV at satTiration density 
Po = 0.16 fm~^, we also use as a constraint the EOS of high density nuclear matter determined 
from analyzing the collective flow of relativistic heavy-ion collisions [45]. Here, we use the same 
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scaling functions of meson masses and meson-nucleon coupling constants as for the model SLl 
constructed in Refs. [29, 30]. The model SLl, however, does not have the vacuum renormaliza- 
tion. The resulting new parameter set RSLl and the corresponding saturation properties are 
tabulated in Table 1. Compared to the parameter set SLl, in the RSLl the coupling constant 
of the oj meson is significantly reduced due to the vacuum contribution. Consistent with the 
finding in Ref. [37], the nucleon effective mass at saturation density increases significantly due 
to the vacuum correction. This means that the in-medium nucleon feels a scalar repulsion from 
the Dirac sea, and it prevents nucleons in the Fermi sea from falling to the fully-occupied Dirac 
sea by sustaining a larger mass gap. 

Table 1: Parameter sets RSLl, RSL2 and SLl. The vacuum hadron masses are M = 938MeV, 
rUfj = 590MeV, m^j = 783MeV and rup = 770MeV. The coupling constants given here are those 
at zero density. The symmetry energy is fitted to 31.6MeV at saturation density. 
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Figure 1: (Color online) The pressure for different models versus density in symmetric nuclear 
matter. The shaded region is given by experimental error bars [45]. 
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Shown in Fig. 1 are the EOSs of symmetric nuclear matter for the RMF model SLl and the 
relativistic Hartree models. The shaded area indicates the experimental constraint on the high- 
density EOS [45]. The model parametrization in the SLl was selected such that the resulting EOS 
passes through the middle of the constraint. Comparing to the EOS with the SLl, it is seen that 
the pressure with the RSLl is much reduced at high densities due to the vacuum correction. In 
order to bring the predicted pressure closer to the experimental region, we use another parameter 
set RSL2 that is adjusted to have a larger nucleon effective mass at saturation density, as listed 
in Table 1. As shown in Fig. 1, the moderate readjustment of parameters results in a significant 
increase of pressure at high densities where the rearrangement term plays an important role. 

We have noticed that the energy density at high densities is dominated by the repulsion 
provided overwhelmingly by the to meson. This is the main origin of the hindrance toward a 
chiral symmetry restoration observed using energy densities predicted by most meson-exchange 
many-body theories. In the simple RMF models, the repulsion increases linearly with density. 
After taking into account the higher-order correlations beyond the mean-field approximation 
within the DB approach, though the EOS at high densities was softened, the increasing tendency 
of the in-medium chiral condensate was not overturned [18, 19]. On the other hand, this may also 
be understandable since these approaches arc not really expected to give a correct description 
of chiral condensate at least near the critical density for the chiral symmetry restoration. It 
was found very recently within the Brueckner approach with chiral limit that nucleon-nucleon 
interactions lead to a weaker reduction of the in-medium chiral condensate in the density region 
where the effective field theory is applicable [26], while little attention was paid to the saturation 
properties. The decreasing tendency of the in-medium chiral condensate can be obtained in RMF 
models by including the density-dependent coupling constants, namely the density-dependent 
ratios and C^} [20]. In this work, instead of modifying the ratios Ca and C^^, we will investigate 
the role played by the vacuum corrections. 

In Fig. 2, we show the in-medium chiral condensate evaluated using Eq.(19) with the RMF 
models NL3 [46] and SLl and the relativistic Hartree models RSLl and RSL2. Compared to the 
NL3 result, the in-medium chiral condensate with the SLl is largely reduced at higher densities. 
This can be attributed to a much softened EOS with the SLl [29], compared to that with the 
NL3 [45]. At high densities, however, the chiral condensate with the SLl increases with density 
given the fact that the constant ratio in SLl does not suppress the vector part of the energy 
density. The latter is quadratic in density, also see Eqs.(9) and (19). It is seen that the mass 
dropping given by the BR scaling is not sufficient to reduce the chiral condensate to zero due to 
the nucleon-nucleon interactions in nuclear matter. 
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Figure 2: (Color online) In-medium chiral condensate for various models. The "vacuum" is 
calculated with the RSL2. 

As the chiral condensate arises from the vacuum, the appropriate description of the Dirac 
vacuum contribution is imperative. The curve in the figure denoted by "vacuum" is calculated 
from the vacuum term Af" in Eq.(9). It is seen that the vacuum contribution to the in- medium 
chiral condensate becomes increasingly more important at higher densities. Moreover, the vac- 
uum correction in the energy density also includes the part that modifies the nucleon effective 
mass. In the renormalized model RSLl the chiral condensate further goes down with the increas- 
ing density. We note that the vacuum energy density /S.8 is positive and the scalar field a that 
provides a big attraction is reduced due to the vacuum correction. Thus, the vacuum correction 
actually gives rise to a scalar repulsion that impairs the role of the vector repulsion provided by 
the uj meson. As a result, the much reduced g^) with a consistent but smaller reduction of g^-, 
subject to the given saturation properties, results in a softened EOS and brings down the chiral 
condensate at high densities in the RSLl, compared to that with the SLl. 

Interestingly, when we use the parameter set RSL2 that is adjusted to have an EOS more 
close to the experimental constraint at high densities, see Fig. 1, the in-medium chiral condensate 
decreases further to zero by coincidence. As shown in Eqs.(9) and (10), the vacuum correction 
Af" contributes oppositely to the energy density and the pressure. This provides the intrinsic 
mechanism for the different changing tendencies in the energy density and the pressure, intrigu- 
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ing the shift going from the RSLl to the RSL2 as observed in Figs. 1 and 2. Numerically, we 
find that a larger cancellation of contributions from the scalar and vector mesons in the RSL2, 
compared to the RSLl, is responsible for the big half of the reduction of the chiral condensate 
from the RSLl to the RSL2, whereas the rest reduction arises directly from the variation of the 
vacuum correction. This further shows a direct link between the in-medium chiral condensate 
and the vacuum correction from the Dirac sea. With the vacuum correction, the nucleon effec- 
tive mass drops more slowly in nuclear matter. As it was also stressed in the modified RMF 
models [20] and the Brueckner approach with the chiral limit [26] , there is a decoupling between 
the reductions of the nucleon effective mass and the chiral condensate at high densities. Actu- 
ally, this decoupling exists in most hadronic approaches where the in-medium chiral condensate 
increases whereas the nucleon effective mass drops at high densities. 

In this work, we have concentrated on investigating effects of the vacuum corrections on the 
in-medium chiral condensate. However, if the chiral condensate does not drop sufficiently at 
high densities, other corrections, such as, high-order correlation contributions, may help further 
bring it down. Another simple way is to reduce the in-medium vector coupling constant. The 
expectation of vanishing chiral condensate at high densities would actually add some constraints 
on the in-medium scaling of the vector coupling constant in the framework of the BR scaling. 
It is also interesting to point out that the parameter used in the scaling function for the meson 
masses is consistent with that extracted from recent low-density experimental data, see Ref . [29] . 
However, it may be risky to apply the currently used scaling function having a linear density 
dependence of meson masses to high densities. In Rcf. [32], different scalings were suggested 
for the regions below and above the saturation density. Moreover, one may consider the use 
of the nucleon mass scaling. Theoretical uncertainty may also exist in the approximations to 
determine the sigma terms for the vector and scalar mesons. Certainly, all these details deserve 
further investigations. Nevertheless, regardless of the very details and the uncertainties that ebb 
and flow, a definite fact is that the vacuum correction brings out a new mechanism favoring the 
reduction of the in-medium chiral condensate. 

5 Summary 

In summary, we have studied the vacuum corrections in the relativistic Hartree model using 
the density-dependent parameterizations according to the BR scaling. The in-medium chiral 
condensate have been investigated using the Hellmann-Feynman theorem together with phe- 
nomenological relations for the current quark mass derivative of hadron masses. Constrained by 
the saturation properties, the vector repulsion is reduced by the scalar repulsion arising from 
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the Dirac sea. This results in the softened EOS and the reduction of the in-medium chiral 
condensate at high densities. Also, the vacuum correction from the Dirac sea plays a direct 
and important role in generating the in-medium chiral condensate. The interplay between the 
vacuum corrections to the energy density and the pressure provides a mechanism favoring the 
chiral symmetry restoration in dense matter. 
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